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FOREWORD 

This  paper,  The  Accuracy  of  Orbit  Prediction  of  Nqr - 
Cooperative  Earth  Satellites,  toy  R.  L.  Dunconbe  of  the  U.  S, 
'Kaval  Observatory >  is  one  o?  a  group  of  papers  prepared  as 
a  result  of  a  meeting  of  OEG  consultants  and  other  experts 
held  on  18  January  1962  to  assist  NAVWAG  in  its  study  of 
certain  aspects  of  the  threat  of  hostile  spacecraft  to  naval 
operations. 

A  resume  of  the  discussions  at  the  meeting  has  already 
been  published- 

Summary  Report  of  the  PEG  Consultants  Conference  on 

the  Threat  of  Hostile  spacecraft  to~Nav al  Operations 

"Secret ' 

and  additional  papers  will  be  published  shortly  dealing  with 
certain  sensors  in  military  satellites  and  with  the  accuracy 
of  prediction  of  the  orbits  of  active  and  passive  satellites. 
No  attempt  has  been  made  to  reconcile  the  contents  of  the 
group  of  papers  as  a  whole. 

This  document  does  not  necessarily  represent  the  opinion 
of  OEG,  NAVWAG,  or  the  U.  S.  Navy. 
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ABSTRACT 

Two  systems  for  tracking  passive  (non-cooperating)  satellites 
are  considered:  the  Navy  Space  Surveillance  System  (NavSpaSur) 
and  the  "Millstone"  active  radar.  Accuracy  estimates  are 
presented  for  both.  Future  improvements  to  both  systems  and 
combinations  of  them  are  also  considered.  The  effect  of 
large  solar  flares  on  prediction  accuracy  is  estimated. 
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THE  ACCURACY  Of  ORBIT  PRFf>;C7?0N  OF 
NON-COOPERATIVE  EARTH  SATEI.ET7ES  sU» 


At  the  present  time  t.v  prirc  ipfe .  jotit; r  of  non-cooperative 
earth  satellites  is  the  U  S  Naval  Space  Sur veillar.ce  System, 
(reference  (a))  an  electronic  sensor  spanning  the  southern 
United  States,  This  system  comprises  several  CW  transmitters 
and  interferometric  receivers,  and  a  central  observation 
analysis  and  orbit  prediction  facility  The  system  is  able 
to  monitor  earth  satellites  having  an  orbital  inclination  of 
approximately  32  or  greater  The  principal  drawback  of  the 
system  is  that  it  is  statically  oriented  with  respect  to  the 
surface  of  the  earth  and  hence  must  depend  on  the  rotation  of 
the  earth  to  carry  it  to  the  orbital  plane  For  the  same  reason 
its  view  is  constrained  to  two  regions  of  orbital  longitude 
which  change  only  because  of  the  secular  variations  of  the  node 
and  perigee  This  restriction  places  definite  limitations  on 
the  "one  pass"  capability  of  the  system 

Of  the  other  present  methods  of  monitoring  non-cooperative 
earth  satellites  probably  the  Millstone-type  tracking  radars 
are  the  most  effective  These  steerabie  radars  are  capable  of 
observing  satellites  over  an  extended  arc  of  orbital  longitude, 
thus  providing  better  "one  pass"  capability  than  the  NavSpaSur 
system.  The  principal  restriction  with  this  radar  however  is 
the  effective  area,  of  search,  which  is  so  limited  as  to  make 
the  cost  of  complete  sky  monitoring  by  such  equipment  prohibitive. 

The  accuracy  of  orbit  prediction  is  limited  by  several  factors 
a)  errors  of  observation  (including  instrumental  errors,  refraction 
error,  timing  error,  station  errors,  poor  distribution  in  orbital 
longitude,  etc  ):  and  b)  inadequacies  of  theory  arising  from 
inability  to  anticipate  variations  in  the  atmospheric  drag 
perturbation,  inaccuracies  of,  or  neglect  of  terms  in  the 
disturbing  potential,  and  uncertainties  in  the  adopted  astro¬ 
nomical  and  geodetic  constants 

Probably  the  largest  body  of  observational  data  ^although 
not  the  most  accurate  due  to  instrumental  limitations*  is  that 
collected  by  the  NavSpaSur  system  Approximately  140  satellites 
are  currently  being  observed  representing  a  cross  section  of 
all  objects  launched  to  date  At  present,  the  elevation  angles 
of  the  NavSpaSur  system  are  recorded  in  aralog  form  and  can  be 
resolved  to  only  1/20  of  a  degree,  equivalent  to  a  triangulation 
accuracy  of  1  2  miles  at  a  height  of  ) 000  miles  Because  of 
this  limitation  in  observational  accuracy  the  orbit  predictions 
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are  based  on  a  restricted  theory  which  incorporates  only  the 
secular  perturbations  and  in  some  cases  the  larger  periodic 
perturbations  plus  terms  in  the  mean  motion  factored  by  the 
time  an^’  the  time  squared  to  represent  the  mean  atmospheric 
drag  perturbation.  Within  these  restrictions  th_*  following 
statistics  illustrate  the  average  prediction  accuracy  obtain-* 
able  with  the  NavSpaSur  system  for  the  present  earth  satellite 
population.  For  a  sample  of  108  satellites,  (reference  lb)), 
the  error  of  prediction  of  a  crossing  of  the  NavSpaSur  "fence" 
one  week  after  the  last  observation  included  iu  the  orbit 
correction  is  shown  in  Table  I: 

TABLE  I 


Pred.  Error 

No.  of  Satellites 

Sum 

Per  Cent 

0s 

20 

20 

19 

Is 

33 

53 

49 

2s 

23 

76 

70 

3r 

9 

85 

79 

4s 

3 

88 

31 

5s 

4 

92 

85 

6-lCs 

8 

100 

93 

1 1-3 0s 

6 

106 

93 

30s 

2 

108 

At  the  average  height  of  this  sample  (350  miles)  an  error 
of  one  second  is  equivalent  to  slightly  less  than  5  miles 
displacement  in  orbit.  The  NavSpaSur  system  therefore  shows 
an  85%  capability  to  represent  the  positions  of  satellites  in 
orbit  one  week  in  advance  within  25  miles.  Due  to  the  low 
average  height  of  the  sample  and  the  fact  that  many  of  the 
objects  (such  as  burned  out  rocket  casings,  nose  cone  shields, 
etc.)  have  a  high  drag  cross-section  to  mass  ratio,  the 
principal  reason  for  this  prediction  error  seems  to  be  the 
inability  to  represent  the  atmospheric  drag  perturbation, 
especially  in  the  last  stages  of  crbital  iecay. 

The  uncertainty  of  predictions  based  or.  observations 
limited  in  orbital  longitude  is  illustrated  by  the  graphs 
In  flguru  l,  where  predictions  based  on  NavSpaSur  observations 
only  are  compared  with  NASA  observations  well  distributed  in 
orbitrl  longitude,  (reference  (b)).  In  both  cases,  however, 
about  90%  of  the  observed  lines  of  sight  agree  with  the 
predictions  within  three  tenths  of  a  degree. 
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FIG.  1:  DISTRIBUTION  OF  ANGULAR  DIFFERENCE  e  BETWEEN  OBSERVED 
LINES  OF  SIGHT  FROM  SAG  STATIONS  AND  PREDICTED  LINES  OF  SIGHT 
FROM  NA VSPASUR  ELEMENTS  FOR  I960  IOTA  2 
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Another  body  of  observational  orbit  data  now  being  analysed 
Is  that  from  the  NASA  Mercury  Project.  Observations  are  made 
with  ITS  16  and  Verlort-type  radars  of  a  transponder  mounted 
aboard  the  Mercury  capsule.  An  analyse  of  the  quality  of  the 
MA-5  real-time  orbit  determination  by  tbe  Mercury  Project 
tracking  and  computing  complex  by  Schiesser,  (reference  (c/) 
indicates  the  following: 

"a.  The  uncertainty  in  position  appears  to  grow  linearily 
with  time  for  the  limited  period  covered  by  MA-5  data. 

"b.  The  FPS  16  is  roughly  twice  as  good  as  the  Verlort- 
type  radar  for  orbit  determinations  based  on  data  from  a 
single  radar. 

"e.  There  is  a  useful  correlation  between  the  precision 
of  single  radar  orbit  determinations  and  the  maximum  elevation 
of  vhe  vehicle  with  respect  to  the  radar,  which  can  be  used  to 
assess  precision  as  a  function  of  tracking  conditions. 

"d.  Single  radar  orbit  determinations  based  on  data  obtained 
to  a  maximum  elevation  much  below  ten  degrees  is  nearly  useless 
for  accurately  predicting  the  future  position  of  a  vehicle, 

"e»  Increasing  the  number  of  participating  radars  beyond 
any  given  number  will  not  always  result  in  greatly  improved 
precision  because  of  the  influence  of  other  important  factors 
For  MA-5  the  greatest  improvement  in  precision  occurred  when 
the  number  of  participating  radars  was  increased  from  two  to 


he  concludes  that  "when  predicting  ahead  using  a  set  of 
orbital  elements  determined  in  real  time,  the  position  of  a 
vehicle  will  be  known  to  lie  within  a  sphere  of  radius  U 
For  predictions  up  to  two  orbits  ahead,  it  has  been  determined 
that  U  may  be  approximately  by  the  expression: 

U  -  A  +  B.t 

where  t  (minutes)  is  the  elapsed  time  from  the  time  of  the 
data  used  in  the  computation  of  the  set  of  orbital  elements 
to  the  time  of  interest," 

For  predictions  based  on  observations  from  one  radar  only 
the  value  of  A  is  approximately  700  yards  for  the  ITS  16  and 
about  1600  yards  for  the  Verlort-type  radar.  For  single  station 
observations,  tracking  to  a  maximum  elevation  of  fifty  degrees, 
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B  ip  found  to  be  nearly  65  yards  per  minuce  for  an  FPS  16  and 
about  130  yards  per  minute  for  a  Verlort-type  radar.  This  leads 
to  radii  of  uncertainty  of  3.7  miles  and  7.6  miles  respectively 
for  predictions  one  orbit  in  advance  of  the  observations  made 
by  the  FPS  16  or  Verlort-type  radars  For  predictions  based  on 
observations  by  three  stations,  the  mean  values  of  A  and  B  are 
650  yards  and  15  yards  per  minute  respectively,  leading  to  a 
radius  of  uncertainty  one  orbit  in  advance  of  1.2  miles. 

These  figures  must  be  considered  to  be  the  maximum  prediction 
precision  attainable  with  FPS  16  radar  for  a  non-cooperative 
satellite  of  similar  size  and  orbit.  It  is  likenv,  due  to  the 
increased  uncertainty  of  observations  made  withou  the  aid  of  a 
transponder,  that  the  actual  prediction  accuracy  will  be  degraded 
by  one  or  two  orders  of  magnitude. 

An  extensive  correlation  of  changes  in  the  geomagnetic 
planetary  index,  as  a  result  of  solar  flares,  with  the  atmospheric 
drag  perturbation  effect  on  seven  different  earth  satellites  has 
been  made  by  Jacchia,  (reference  (d)).  The  resulting  accelerations 
are  illustrated  in  Figure  2,  taken  from  that  report.  This 
illustrates  that  for  a  satellite  of  similar  mass  and  drag  cross- 
section  to  Vanguard  1,  at  a  height  of  400  miles,  an  orbit 
prediction  fitted  to  observations  in  the  period  Nov,  12.5-13,5 
and  having  a  mean  epoch  of  Nov.  13.0  would  be  found  not  to 
represent  an  observation  at  Nov,  14.5  by  slightly  over  two-tenths 
of  a  mile  from  this  cause  alone.  It  should  be  noted  that  the 
time  between  onset  of  the  flares,  as  indicated  by  the  arrows  in 
tt  time  argument,  and  the  change  in  dP  or  the  geomagnetic 

dt 

planetary  index,  is  variable  and  presently  unpredictable. 

The  number  of  solar  flares  of  similar  intensity  (equal  to 
or  greater  than  2+)  recorded  since  mid-1958  are  shown  in  Table  II, 
The  mean  sunspot  cycle,  with  which  the  incidence  of  these  flares 
is  correlated  is  shown  in  Figure  3. 
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FIG.  2:  ATMOSPHERIC  DRAG  OF  SEVEN  ARTIFICIAL  SATELLITES  DURING 
THE  NOVEMBER  1960  EVENTS,  COMPARED  WITH  THE  GEOMAGNETIC 
PLANETARY  INDEX  ap  AND  THE  SOLAR  FLUX  AT  10.  7  cm. 

FOR  DETAILED  EXPLANATION,  SEE  TEXT. 
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TABLE  II 


Solar 

Flares  of 

Impor  tance 

2  * 

' 

1955 

1959 

1960 

1961 

Jan. 

14 

2 

- 

Feb. 

17 

8 

- 

Mar. 

14 

6 

1 

Apr. 

20 

11 

2 

May 

24 

9 

1 

June 

17 

16 

5 

July 

31 

28 

6 

24 

Aug. 

29 

15 

13 

2 

Sept, 

22 

8 

6 

2 

Oct . 

20 

- 

a 

V 

- 

Nov, 

9 

12 

14 

- 

Dec. 

15 

4 

1 

1 

126 

173 

101 

38 

The  foregoing  evidence  has  been  presented  to  indicate  the 
accuracy  with  which  the  orbit  of  a  non-cooperative  satellite 
sight  be  predicted,  using  first,  NavSpaSur  observations 
extending  over  several  days  and  second,  single  pass  observations 
of  NASA  Mercury  Project  collected  frcra  one  or  three  stations. 

Improvements  in  instrumentation  and  theory  will  undoubtedly 
occur  in  the  next  few  years  to  greatly  increase  the  ca~ ability 
of  detecting  and  predicting  the  orbits  of  non-cooperative 
satellites.  Kozai,  (reference  (e)),  has  summarized  the  most 
recent  values  of  the  constants  associated  with  the  zonal  and 
tesseral  harmonics  in  the  potential  of  the  earth  as  derived 
from  the  orbital  variations  of  earth  satellites.  Further 
improvement  in  these  and  the  geodetic  constants  may  be  expected 
from  Project  Anna,  Probably  the  largest  group  of  observations 
from  which  improved  constants  may  be  derived  are  those  collected 
by  Project  Trp.nsit,  altnojgh  the  numerous  photographic  observa¬ 
tions  of  the  Smithsonian  Astrophyslcal  Observatory  tracking 
stations  will  be  of  value  when  precisely  reduced.  It  may  be 
noted  that  the  next  few  years  offer  an  excellent  opportunity 
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for  collection  of  precise  observations  with  a  minimum  of  solar 
flare  interference  (see  Figure  3),  The  most  recently  determined 
values  of  astronomical  and  geodetic  constants  have  hesn  reviewed 
hy  Kaula  (reference  (f))  It  is  important  to  the  future 
combination  of  observations  and  derivation  of  improved  constants, 
tb  vt  this  or  come  other  mutually  agreeahle  set  of  constants  be 
adopted  for  use  hy  all  earth  satellite  ohserving  facilities, 


Refinement  of  the  presently  used  model  of  the  earth's 
atmosphere  must  await  the  collection  of  earth  satellite 
observations  at  higher  altitudes  and  for  various  drag  cross- 
section  to  mass  ratios.  Application  of  presently  observed 
atmospheric  drag  effects  to  satellites  outside  the  region  of 
observation  or  of  differing  configuration  may  he  subject  to 
error  because  of  presently  unknown  variations  of  atmospheric 
parameters.  The  atmospheric  drag  effect  and  unpredictable 
variations  thereof  will  continue  to  be  important  to  the  accuracy 
of  long,  and  short,  range  predictions  of  non-cooperative 
satellites;  particularly  those  resulting  from  enemy  action, 
which  will  he  of  greatest  concern.  It  seems  impossible  to 
anticipate  a*  what  height  to  expect  non-cooperative  enemy 
satellites  While  higher  altitudes  would  promote  orbital 
stahility  for  surveillance  purposes,  the  trajectory  of  a 
missile  from  such  a  high  satellite  hack  to  the  surface  of  the 
earth  would  he  long  enough  to  allow  detection  and  defensive 
action.  On  the  other  hand,  for  reasons  of  booster  limitations 
or  short  mission  life-time,  non-cooperative  enemy  satellites 
might  he  placed  in  relatively  low  orhits.  Under  any  circum¬ 
stances,  no  precise  Information  as  to  the  mass  or  drag  cross- 
section  of  enemy  satellites  will  be  availahle  and  hence 
prediction  of  their  orhits  will  always  be  affected  hy  the 
inaccuracy  of  our  approximations  to  these  parameters. 

Current  improvements  in  the  power  and  instfumemtation  of 
the  NavSpaSur  system  will  allow  detection  of  targets  of  one 
square  meter  at  a  height  of  2500  miles  with  a  triangulation 
accuracy  of  ,06  miles.  Operation  of  this  detection  and  spotting 
capability  to  give  immediate  alerts  to  steerable  radar  of 
similar  accmccy  in  order  to  obtain  observations  over  long 
arcs  would  seeu  to  provide  the  best  team  for  the  detection  and 
tracking  of  non-cooperative  satellites. 


Incorporating  refined  values  of  the  astronomical  and  geodetic 

constants,  improved  values  of  the  refraction  and  an  extended 

model  of  the  atmosphere  with  a  complete  gravitational  theory, 

should  produce  real-time  orhital  predictions  for  non-cooperative 

earth  satellites  hased  on  these  observations  having  an  accuracy 

one  orbit  in  advance  one  order  of  magnitude  hetter  than  that 

now  attained  with  the  Mercury  tracking  computing  complex, 

« 
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